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The quantum mechanical equivalent of parametric resonance is studied. A simple model of a 
periodically kicked harmonic oscillator is introduced which can be solved exactly. Classically stable 
and unstable regions in parameter space are shown to correspond to Floquet operators with qualita- 
tively different properties. Their eigenfunctions, which are calculated exactly, exhibit a transition: 
for parameter values with classically stable solutions the eigenstates are normalizable while they can- 
not be normalized for parameter values with classically instable solutions. Similarly, the spectrum 
. of quasi energies undergoes a specific transition. These observations remain valid qualitatively for 

' arbitrary linear systems exhibiting classically parametric resonance such as the paradigm example 

^ , of a frequency modulated pendulum described by Mathieu's equation. 
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I. INTRODUCTION 

>...... 

QQ , The highly complicated behaviour of classically chaotic systems does not translate in an obvious and straightforward 

' way into properties of their quantum counterparts. Nevertheless, various features such as energy level statistics and 
\ the spatial structures of wave functions have been identified as more or less faithful quantum indicators of chaos in 
^5 . the classical limit. The extreme sensitivity of classically chaotic systems to the variation of initial conditions appears, 
for fundamental reasons, to figure less prominent in quantum systems. 

Parametric resonance is a useful concept in order to understand the transition from regular to chaotic motion in 
classical dynamic systems. Each of the phase space filling tori of an integrable system (with two degrees of freedom) is 
characterized by a frequency ratio, and the value of this parameter determines the fate of a torus under a perturbation 
according to the KAM theorem Q . When the strength of the perturbation increases, initially small areas of instability 
^ _ proliferate in the form of Arnold tongues. They echo the intricate structure of phase-space regions where motion on tori 
^ and irregular trajectories coexist. Apart from its conceptual importance, parametric resonance has found applications 
^ in systems with many degrees of freedom providing, for example, an important mechanism in the formation of patterns 

FT ™- 

The purpose of this paper is to study the phenomenon of parametric resonance in a quantum context. Hamiltonians 
quadratic in position and momentum are known to provide relevant examples of classical parametric resonance. Since 
for such systems, classical and quantum mechanics of systems are intimately related, the program here is to make 
] explicit what the quantum manifestations of the classical instability look like. In particular, the focus will be on the 
- - ■ quasi-energy spectrum of the quantum mechanical evolution operator or Floquet operator as well on its eigenstates. 
It will be shown that the properties of the quantum system unambiguously reflect the stability and instability of its 
classical counterpart. 

In Section classical parametric resonance is discussed from a general point of view. An analytically solvable 

'A harmonic oscillator 
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model which exhibits parametric resonance on the classical level, is introduced in Section HI: 
is subjected to a perturbation which periodically dilates and squeezes volumes in phase space. The parameter ranges 
for the instability of the classical system are determined. In Section Floquet eigenstates and quasi energies of the 
associated quantum system are calculated in the unstable case. By means of an effective Hamiltonian a comprehensive 
point of view for the possible scenarios is developped in Section ^ Subsequently, the result is shown to persist for 
linear systems with arbitrary periodic frequency modulation. The generalization includes systems described classically 
by Mathieu's equation. Finally, the results are summarized and links to related models are established. 

II. CLASSICAL PARAMETRIC RESONANCE 

Parametric resonance occurs if an appropriate parameter of a classical dynamical system is varied periodically in 
time. Stable fixed points of the flow in phase space become unstable for specific values of the period of the parameter 
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variation. As an example, consider a pendulum with mass m and length I under the influence of gravity the support 
of which is moved up and down by an amount Al < I with frequency loq. The equilibrum position of the undriven 
system corresponds to the bob resting vertically below the support. For small oscillations, the pendulum is described 
by Mathieu's equation 

f^+coHt)x = 0, (1) 

with 

Lu"^ {t)=l- Al cos{LUot) > , (2) 

where x{t) is the vertical coordinate of the bob. This equation provides the paradigm for parametric resonance 
in classical mechanics: the behaviour of its solutions has been investigated in great detail |^ as a function of the 
parameters I, Al, and ujq (see Q for an interesting discussion). The equilibrium position is destabilized ioi appropriately 
chosen values of Al and loq- arbitrarily small deviations from the fixed point will grow exponentially in time. This 
is characteristic for parametric resonance. For other parameter values the elliptic fixed point of the undriven system 
survives the perturbation. Smooth boundaries decompose the parameter space into regions with stable and unstable 
behaviour. On the separating boundary, a third marginal type of behaviour occurs which interpolates between stable 
and unstable motion. 

Parametric excitation can also have the opposite effect which is less intuitive: an unstable fixed point may become 
a stable one. Indeed, the inverted pendulum turns stable under specific conditions on and ^{t) a phenomenon 
known as parametric stabilization. Small deviations from the vertical position do not grow without bound but the 
pendulum remains in the neighbourhood of its initial position for all times. 

Qualitatively, these phenomena show also up for other types of periodic driving, + T) = (^'^{t), with a modified 
division of parameter space into stable and unstable regions. In the following, a particularly simple quantum system, 
known to exhibit parametric resonance classically, will be studied in detail. 



III. A SIMPLE MODEL 



Consider a harmonic oscillator with frequency uj, Hq = /2m + muj'^x^ /2, subjected to a periodic perturbation 

H{t)^Ho + Hkit), (3) 

where 

Hk{t) = ^{xp + px)ST{t) , a € IR , (4) 
with an infinite comb of delta functions, 

oo 

5T{t)= S{t~nT). (5) 

In between the times r„ = nT, n ^7L, the particle with mass m moves freely in a quadratic potential, while at times r„ 
it experiences an impulsive kick with amplitude axp. In fact, the Hamiltonian does not depend on three but only two 
parameters {ljT and a) as is seen from introducing r = t/T, and rescaling simultaneously p — *■ ^Jmlop, q — > q/ ^/nuo. 

From a general point of view, the momentum dependence of the kick is a particular feature of this system. Other 
driven systems which have been studied as models for chaotic motion in both classical and quantum mechanics, such 
as the kicked rotator |9[|, have a position dependent amplitude only. Various authors pointed out the consequences 



which result from this difference 11|, and explicitly solvable models have been studied [ |l2[ . 

Consider the evolution of the system over one period, from t = —e to t = T — e, say, for e ^T. The equations of 
motion are 



d ( x \ _ ( aSxit) 1/m 



dt \ P J \ -muj^ -aSrit) J \ P 



X 



(6) 



The integration over the time interval — e < t < e, is carried out by first introducing a smooth approximation Sj^{t) 
of the delta functions in (p). Then the solution of (0) is given by 
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/ x{+e) \ ^ 



dtM^[t) 



x{-e) 



(7) 



where the symbol T denotes time ordering, and M^{t) is the smooth 2x2 matrix from (^. A straighforward calculation 
gives a simple result in the limit e — *■ 0, 



Mh = lim exp 



+e 



dtM^(t) 



exp 



a 

-a 



e" 
e-" 



(8) 



Consequently, the interaction Hi,{t) rescales instantaneously position and momentum by the amounts exp[±Q;], re- 
spectively. In other words, it generates periodically dilations [Q, which preserve volume in phase space. 

Next, Eqs. are integrated over the remaining time interval between the kicks, e < t < T ^ e. In the limit e ^ 0, 
one obtains, 



Mn = lim 



lim exp 



1/m 
rncj^ 



(T-e) 



cosujT (mcj) "'^sintjT 
—muj sin LoT cosuT 



(9) 



Writing z = {x^p)^ for points in phase space, the time evolution of the kicked harmoic oscillator over one full 
period from t — io t — T~: z{T~) = Mz(0~), is thus due to a kick followed by free motion, generated by the 
matrix 



M = MoMk 



e" cos{ujT) (e"muj) ^ smtoT 
-e"mwsinwT e^^cosojT 



(10) 



being symplectic and of unit determinant. 

Dilations and oscillatory motion compete with each other. The overall character of the motion over one period 
depends on the actual values of the parameters. The matrix M does not vary over phase space. Hence, for given 
values of the parameters w, a and T, all initial values z (7^ 0) evolve in a similar way. The matrix M in (|l|) can 
generate three qualitatively different types of motion, conveniently characterized by their eigenvalues 



A± = -Tr M ± 
2 



-(Tr MY 



cosh a cos ijjT ± V cosh^ a cos^ uiT — 1 



(11) 



The eigenvalues may come as 



(z) a complex conjugate pair with modulus one: \± — exp[±if2]; 
{ii) a real reciprocal pair: |A±| = exp[±/i]; 
{m) a real degenerate pair with modulus one: A_|_ = A_ = ±1. 

If the eigenvalues of M are purely complex [i) , the images of a point of phase space z under the composed action of 
Mo and M^ are all located on an ellipse. Appropriate rescaling of the axes transforms it into a circle such that during 
one period the angular coordinate of a point is seen to increase by an angle u) defined through coscD = coshacoscjT. 
The perturbation changes the original motion only quantitatively ^ i. e. the frequency of the oscillator now is ui instead 

of Ld. 

In case {11) , the real eigenvalues of M indicate a "hyperbolic rotation" parametrized conveniently by the positive 
real number with cosh/x = coshacoscijT. Now the time evolution of the oscillator is changed qualitatively by the 
perturbation Hk{t), from stable to unstable motion. According to the sign of the eigenvalues A± the iterates of a 
phase-space point are located on one (ordinary hyperbolic) or two (hyperbolic with reflection) branches of a hyperbola. 
The parameter regions associated with cases («) and (zz) exhaust the parameter space almost completely. 

Stable and unstable regions are separated by boundaries defined by the condition | coshacoswTl = 1. For the 
parameter values of case (tn), iterates of phase-space points are either on one (A+ = A_ = -|-1) or on two (A+ = A_ = 
— 1) straight lines. This situation may be thought to interpolate between the elliptic and hyperbolic case. 

More precisely, the three cases are related with specific invariant curves in phase space. To see this, associate a 
quadratic form with the matrix M by 
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Q{z) = z^AMz , A = J j . (12) 

In fact, Q{z) is invariant under time evolution over one period: 

Q{Mz) = Q{z) , (13) 

as follows from the symplecticity of M: M'^AM = A. Using Mq and Mk as examples for (i) and {n), respectively, 
one obtains 

Qo(2) oc i (p^ + m^tj^x^) , Qk{z) oc px oc {p + x)'^ - {p - x)'^ , (14) 
while the marginal situation (in) implies 

Qtiziz) cc {p±x)'^ , (15) 
using coswT = ±1/ cosha and sinwT — tanha. These results can be summarized by writing 

Q{z)^^{P^ + n^X') , (16) 

where {X,P) are related to {x,p) by appropriate linear canonical transformations. Then, a variation of the original 
parameters is reflected in a change of the factor il? multiplying the quadratic term in X: f2 is either real, purely 
imaginary, or zero. In physical terms, the time evolution of the system over one period is effectively that of a particle 
in an attractive or repulsive quadratic potential, or in the absence of a potential. Not surprisingly, the quadratic 
forms Q{z) will play an important role for quantum parametric resonance. 

If coso-iT = 1, or ujT — 27rfc, fc e 2, an arbitrarily small kick amplitude a renders the system unstable. This 
particularly simple situation — the harmonic evolution hasno net effect: Mq = 1, — is called resonant. Similarly, for 
ujT — 27r(fc + 1/2), fc e 2, the time evolution is resonant with reflection: each point z in phase space is mapped to 
(— z) by the harmonic time evolution: Mq = —1. 



IV. QUANTUM PARAMETRIC RESONANCE 



The quantum mechanical harmonic oscillator with an impulsive force 



^ ^ 2m 



x^ 



-{xp + px)5T{t) 



is described by the Hamiltonian operator 

(17) 



Position and momentum operators x and p satisfy the fundamental commutation relation [x,p\ — ih. The long-time 
behaviour of the system is determined by the Floquet operator = U(tQ + T^to) = Ut- It maps a state at time tQ 
over one period to another state: J^|V'(io)) = |V'(*o + T)). Choose the time io = Tq — 0^ just before the kick at 





t = 0. Then, due to the ^-type interaction the Floquet operator is a product of two unitary operators 

^ = UqUu . 



The first unitary generates harmonic motion with frequency lo from t 

i-T~ 



C/o(r-,0+) =r / dtexp -iH{t)/h 



0+ 



exp 



0+ to i = T- 



P 

2m 



T 



(18) 



(19) 



where T denotes time ordering again. The second unitary describes the effect of the kick at time t ^ 0. As before, 
the time-ordered product is evaluated explicitly by approximating the S distribution as a strongly peaked, smooth 
function Sfp and taking the limit e ^ |12], 



C/fc(0+,0-) =T 



dt exp 



'iH{t)/n 



exp 



2% 



[xp + px) 



(20) 



As shown in |15|, the resulting "squeeze operator" |14| has an implementation in a quantum optical context. Some 
properties of the operator in the exponential in ( pO| ) have been studied in |l6[ . Using the shorthand z = (x,p)"^, the 
action of the Floquet operator T on position and momentum operators works out as in the classical case: 
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z{T-) = Tz{0-)T'' = MoMkz{0-) = Mz(O-) , 



(21) 



where 

z{T-) = Uoz{0+)U^o = Moz{{)+) , (22) 
z{{)+) = Ukz{Q-)Ul = Mkz{i)-) . (23) 

These relations fohow from either integrating the hnear Heisenberg equations of motion or from expanding and 
resumming the exponentials involved. On comparing Eqs. (^,^) to the result (|2^), the intimate relation between the 
classical and quantum time evolutions generated by a quadratic Hamiltonian is evident. 

Therefore, the quantum system is expected to inherit the division of parameter space into stable and unstable 
regions on the basis of Eqs. (|2^), i.e. through the properties of the matrix M. How do the different parameter 
regions manifest themselves in the quantal framework? When moving in parameter space from a classically stable to 
an unstable region, the eigenstates |^^) as well as the quasienergies of the Floquet operator, 

=cxp(-*i?^T/?i)|0^), (24) 

will undergo qualitative changes. The merit of the present model is that one can determine and explicitly 
for characteristic cases. This will be done in Section 5 by means of an effective Hamiltonian _ffgfj which generates the 
time evolution over one period. 

Before presenting the unified treatment, the most interesting case of a Floquet operator associated with a classically 
unstable region of parameters will be studied in detail. Consider, for simplicity, the resonant case, defined by ljT = 
27rfc, fc € Z: the classical motion becomes unstable for any nonzero value of a. The Floquet operator T reduces to 
just the kick operator Uk. It is the simplicity of the action of the kick on position and momentum eigenstates which 
allows one to construct the complete set of ((5-normalized) eigenstates of the time evolution operator. 

According to ([2^), the position operators x± just before and after the kick at time t — Q are related by 

x+ = UkX-Ul^e°'x_, (25) 

and the eigenvalue equations read 

x+\x+) — x+\x+) , x^\x^) — x-\x-) . (26) 
Multiplication of the second equation with Uk from the left and using (^5|) implies that 

x_ {Uu\x^)) = e-"x^ {Uk\x^)) . (27) 

Hence, the state Uk\x-) is an eigenstate of in other words, the operator Uk maps a position eigenstate |a;) to 
another position eigenstate with eigenvalue e~°'x: 

Uk\x) ^ ^=\e-^x) . (28) 



The factor l/\/ p{x) accounts for a possible change of the normalization of the states. It guarantees the completeness 
relation to hold in the form dx\x){x\ = 1 (cf. [0), and it is fixed by the following condition: 

5{x - x') = {x\x') = {x\ulUk\x') = -^(e-"a;|e-"a;') 

pyx) 

^ Sie-^ix-x'))^ ^5{x-x'). (29) 



p{x) p{x) 
This requires p{x) — e" leading to 

C/fc|x) =e-"/2|e-"x), (30) 
which agrees with the result in [ p^ . The action of U^ on a momentum eigenfunction \p) reads 

C/fcb) =e"/2|e», (31) 
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as follows from a similar argument for p or from exploiting {p\x) = {p\UlUk\x) , etc. 

The transformation property ( ^o|) provides the clue to determine the Floquet eigenstates as linear combinations of 
position eigenstates. Imagine to iterate a position eigenstate |a;o) forwards and backwards by applying powers of Uk 
and U^^ = t^fcj respectively. After multiplying with appropriate phase factors candidates for eigenstates oi T = Uk 
have the form: 



xo,^i) = ^ V e*^"e-""/^|e-""xo) . (32) 



1 —an 



n— — oo 



As it stands, \x{)^^) is indeed an (improper) eigenstate of !F with eigenvalue exp[— z/i] since 

'271 

n— — QO 

1 



2tt ^ 

n— — OO 



= e-'^\xf,,fi). (33) 

There is no restriction on the values of /i G [0, 27r)-hence the spectrum of quasi energies is continuous. As one can 
check explicitly, the states \xo,fj) do provide an orthonormal complete set in Hilbert space if the label xq is from 
either one of the two intervals [1, e") or [— e", —1), a > 0, say, 

oo 

{xo,fi\x'„^i') = — Yl e-^''"e''^'-e-"("+™)/2(:Eo|((7l)"t/rko) 

m.n— — OO 
n,m— — OO 

OO 

= — y e^(^'-^)"e-""(5(e-""(xo-x;,)) 
27r — ' 

n— — OO 

= 5{^i - ^i')S{xo - x'o) . (34) 

Note that for xq and Xg chosen from the same the reference interval, the scalar product of position iterates with 
different n and m vanishes and, thus, gives rise to a Kronecker delta 6nm- 

Contrary to one's first expectation, the continuous label Xq € [1, e"] does not imply a continuous degeneracy of the 
quasi energies Ej but only a countable one. This is due to the fact that eigenstates of the position operator x are 
improper ones, i.e. not normalizable. To see this in a simple way, consider a particle in a one-dimensional box. Its 
energy eigenstates span a Hilbert space with countable dimension. Nevertheless, one needs continuously many states 
to span this space by (improper) states \x). In Section 5.3 the countable degeneracy of the quasi energies will also 
follow from general considerations. 

The completeness relation reads 

d,Jdx\x,,){x,,\+Jd, J dx\x,,){x,,\ + lu^Jdx\x)ix\^l (35) 

1 0-1 -e 

where the integration extends over both fundamental intervals, and the state \x = 0) has to be included as an 
additional eigenstate of the Floquet operator. This state is clearly orthogonal to all other eigenstates of T in Eq. 

The Floquet operator associated with a classically instable region is thus seen to have a continuous spectrum of 
quasi energies, and the supprot of its eigenfunctions in configuration space is not compact but extends to infinity. In 
addition, there is one state localized at the origin, x = 0, thus "sitting on top" of the instable fix point of the classical 
map. 



6 



V. EFFECTIVE HAMILTONIAN 



According to Eq. the Floquet operator J- is a product of two noncommuting operators Uq and Uk which are 
quadratic functions of position x and momentum p. Since the Floquet operator is unitary, on can express it in the 
form = eKp[—iHggT/h], with an "effective" Hamihonian -ffcff- This operator is obtained from entangUng the two 
unitaries into a single exponential using Baker-Campbell-Hausdorff technology as presented in p!9[ | , 



exp 



P 

2m 



T 



exp 



'2h 



(xp + px) 







= exp 





(36) 



In the present case, it is possible to determine explicitly the effective Hamiltonian since the operators in the exponents 
of the product in (pfl) constitute a closed Lie algebra. 



[a;^,p^] =2ifi{xp + px) 
[{xp + px) , p^] = 2ifip'^ , 
[{xp + px) , x'^'\ — ~2ihx^ . 



(37) 
(38) 
(39) 



, equivalent to the "split three-dimensional algebra" |T^. Physical realizations of this algebra can be given in terms 
of creation and annihilation operators of a two-dimensional harmonic oscillator or angular momentum operators. 

Due to the commutation relations, the effective Hamiltonian is necessarily a linear combination of the three quadratic 
generators, 



-HesT — ap^ + hSp' 



-{xp + px) . 



(40) 



The transformation of the operators p and x over one period T through T has been determined in Eq. (pT|), while 
in pod their transformation has been calculated starting from a given quadratic expression (^). Combining these 
results, it is straightforward to determine the effective Hamiltonian: the parameters a, &, c must be expressed in terms 
of the elements Mjk {j, k = 1,2) of the matrix M. After some algebra one finds 



c4f(M.- 



22) 



The real number A is given by 



A 



arcsinhD 
D 



A{D^) e IR 



being thus an even function of the difference between the eigenvalues of the matrix M, 



D = i(A+-A_) = 



1 



-(Afn+M22)'-1 



1/2 



Using ( |l0|) for the matrix elements Aljk, one obtains the expression 



sino^T 



ujT V 2me" 



p- 



-x'^ + uj sinh a cot ujT{xp + px) 



(41) 
(42) 
(43) 

(44) 
(45) 
(46) 



For a = or cj = the relations above reproduce correctly the operators Uq and Uk, respectively. It is interesting to 
note that this result is closely related to the quantized quadratic form, i.e. the expression Q — Q{z) obtained from 
Q{z) in (12) by replacing z ^ z and appropriately symmetrizing, 



Q{z) = fkM'z = cri?eff + tI , 



(47) 



with two constants tr and t. Retrospectively, it is natural to find that the generator for time displacement over one 
period is a function of the quantized quadratic form invariant under the corresponding classical motion. In other 
words, the Floquet operator T commutes with the corresponding quadratic form. 
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Therefore, three types of effective quantum Hamiltonians will show up which are in correspondence with the 
classically stable, unstable, and marginal case. Each of them is unitarily equivalent to one of the following expressions: 

i?cfr « ^ + "'^') ' (48) 

where the operators P and X are unitarily equivalent to x and p, and there are three possible values for the frequency 
51, nonzero real, purely complex or equal to zero. For simplicity, the discussion to follow will take ( p8| ) as a starting 
point. There is no need to use the explicit expressions of H^^ff in terms of the original variables. 



A. Stable elliptic case 

For parameter values associated with a classically stable region, the effective Hamiltonian operator H^q is unitarily 
equivalent to a harmonic oscillator with unit mass and real non-zero frequency il. Consequently, the eigenstates of 
the Floquet operator satisfy: 

T\n) = e-''-\n) , e„ =E„T/h, (49) 
with the familiar oscillator eigenstates \n), 

Hes\n) ^ Er,\n) , E,, = hn (n + , n £ No , (50) 

which provide an orthonormal basis in Hilbert space. The spectrum of quasi energies Sn associated with the operator 
^in(P), 

e„ = fir + mod 27r , n = 0,l,2,..., (51) 

is obtained by 'projecting' the oscillator spectrum En onto the interval [0, 27r). Two possibilities arise: 

1. If or is an irrational multiple of 2tt then all quasi energies Sn are different, e„ ^ Sn'. The union of all Sn is 
dense in the interval from to 27r. The set of eigenvalues {Sn} is countable, and no quasi energy is degenerate. 

2. If QT is a rational multiple of 27r, say r/s, then Sn+s = En for all n, and there are only s different quasi energies: 
eo,£i, • ■ • ,£5-1 e [0, 27r). Each value is countably often degenerate, and the Hilbertspace H decomposes in a 
direct sum. 



H = 0H., (52) 



i/=0 



where the finite-dimensional degenerate subspace Hi, is spanned by the states {{v + ks), k S Mo}- 



B. Marginal case 

Consider now parameter values which imply that the effective Hamiltonian operator is unitarily equivalent to that 
one of a free particle: -ffeff , i.e. il = 0. Its eigenstates are identical to those of a momentum-type operator, 

P\P)=P\P), Pe(-oo,cx)), (53) 

hence require a ^-normalization. Since these states are also eigenstates of the Floquet operator J-', it has a continuous 
set of (5- normalized eigenstates. In the position representation the expansion coefhcients of the states \P) have modulus 
one throughout, and the states do not decay for X —^ ±00. 

The quasi energies eo.fc± = P'^/{2h) take all values in the interval [0,27r). Each quasi energy is countably often 
degenerate: given a state \Pq), all states |Po,fc±) with 



P^, = ±^PS + —2nk, fc=l,2,..., (54) 

are eigenstates with the same eigenvalue: exp[— i iPo,k)^T/{2h)] = exp[-iP§T/{2h)]. In analogy with (||), the Hilbert 
space can be thought of as a continuous product of subspaces each with countable dimension. 
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C. Unstable hyperbolic case 



The quantum system associated with a classically instable region is unitarily equivalent to a particle in an inverted 
quadratic potential, that is, the Hamiltonian in ( ^8| ) with < 0. The potential takes arbitrarily large negative values 
for X ±00, hence the solutions oscillate ever faster for increasing X. No normalizable eigenfunctions exist but 
there are two solutions of Schrodinger's equation for every value of the energy E. The spectrum of quasi energies 
exhibits thus the same features as in the marginal case: the numbers Ef^ take any value between and 27r, and 
each value is countably often degenerate. This is due ot the exponential function which "wraps" the real variable 
E G (—00, 00) around a circle. For the resonant case, the solutions have been given explicitly in Section 4 as a sum 
of delta functions in the position representation with nonuniform amplitudes. 



VI. ARBITRARY FREQUENCY MODULATION 



The results obtained so far have been derived for a model with a particularly simple time evolution. Classically, 
however, the phenomenon of parametric resonance is known to occur for much more general periodic frequency 
modulations. In fact, the behaviour of the associated quantum systems can be shown to exhibit qualitatively the 
same behaviour. Consider the classical Hamiltonian 



muP{t) 2 



(55) 



which generates the equations of motion (|l|). The system, a pendulum with a harmonically modulated suspension 
point, exhibits parametric resonance as shown by Mathieu for the choice of a;^(i) in (^). Write the Floquet operator 
of the system as a product of N unitary operators for n consecutive time intervals of length T/N, 



J- = Texp 



dtH{t) 



N-l 



Yl U{tn+l,tn), 



tN=T- 



In the limit N s- 00 the length of the time intervals goes to zero, and one has approximately 



Uitn+l.tn) = Texp 



exp 



dt 



2m 



2m 



(56) 



(57) 



T 
N 



where the number takes a value between u;^(i„+i) and w^(i„). Being quadratic in position and momentum, any 
two adjacent exponentials can be entangled by means of a Baker-Campbell-Hausdorff relation. The result is another 
exponential bilinear in x and p as follows from the algebraic properties (^). Repeating this process for ever larger 
values of N, the Floquet operator tends to 



T — exp 



— yup + vx + — [xp + px 



= exp 



(58) 



According to the values of the parameters u, w, and w, the effective Hamiltonian necessarily will be (unitarily 
equivalent to) one of the three possible types discussed exhaustively in Section 5. As long as the frequency modulation 
is periodic, no qualitatively different behaviour will occur within this class of systems. Clearly, the separation of the 
parameter space into stable and unstable regions will depend subtly on the actual function Lo'^{t). 



VII. DISCUSSION AND OUTLOOK 



The main result of the present paper is the insight that the concept of classical parametric resonance has a well- 
defined quantum analogon. Classically, in a periodically driven linear system three qualitatively different types of 
motion are possible, according to the chosen parameter values. This division of parameter space is mirrored quantum 
mechanically by qualitatively different spectra of the Floquet operator. Its eigenfunctions associated with a classically 
stable region are normalizable, and they become singular when moving into the classically instable parameter region. 
This has been made explicit by constructing them in terms of (improper) position eigenstates. 



9 



The class of time-dependent harmonic oscillators is well-known to provide insight into various aspects of classical 
and quantum mechanics. Therefore, some of the results can be found implicitly in work dealing with other properties 
of such systems. Closest to the present approach is, maybe, work by Perelomov in a discussion of coherent 

states, a group-theoretical approach to driven linear systems is presented based on work in [^iy22[ . Further, the paper 
I p^ deals with parametric resonance and its generation within the realm of quantum optics. In this context, is 
also interesting, where exponential divergence of the energy expectation value has been derived for an oscillator with 
periodically modulated mass — however, no global picture of quantum parametric resonance has been established. 

The impact of an instaneous mass change of a quantum harmonic oscillator also gives rise to squeezing, and the 
impact on the variance of position, for example, has been studied in p4[ . Non-periodic frequency modulations have 
been studied by various authors, starting with the now well-known Caldirola-Kanai oscillator with exponentially 
increasing mass which, effectively, gives rise to damped motion of the oscillator. Among other work, exact 

solutions have been found for a polynomial time-dependence, Lo'^(t) — uP't^, 6 > in pTj . 

Finally, it is worthwhile to point out that the periodically driven oscillator systems studied in the present paper 
suggest a natural generalization: modify the amplitude of the kick according to 

H{t) = ^0 + I {V{x)p + pV{x)) Srit) , (59) 

where V{x) is some smooth function. Such a kick can be shown to generate nonlinear maps of configuration space 
onto itself It will be shown elsewhere that features such as quasiperiodicity and a devil's staircase, which are 
characteristic of classically chaotic dynamics, also play a role in the time evolution of such a quantum system. 
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